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FOREWORD 


This  report  was  prepared  in  the  Mechanics  and  Surface 
Interactions  Branch  (AFWAL/MLBM) ,  Nonmetallic  Materials  Division, 
Materials  Laboratory,  Air  Force  Wright  Aeronautical  Laboratories, 
Wright-Patterson  AFB ,  Ohio.  The  work  was  performed  under  the 
support  Project  No.  2419,  "Nonmetallic  Structural  Materials", 

Task  No.  241903,  "Composite  Materials  and  Mechanics  Technology". 

The  time  period  covered  by  this  effort  was  from  January  to  December 
1980.  Stephen  W.  Tsai  (AFWAL/MLBM)  was  the  laboratory  project 
engineer.  Rodolfo  Aoki  was  a  visiting  scientist  with  MLBM  from  the 
German  Aerospace  Research  Establishment  ( DFVLR) . 

This  is  a  revised  edition  of  the  technical  report  bearing 
the  same  number  published  in  April  1979.  The  asymmetric  laminate 
tapes  are  added  in  this  revision.  The  equations  and  table  numbers 
which  appear  in  the  flow  charts  are  the  same  as  in  Introduction 
to  Composite  Materials,  coauthored  by  S.  W.  Tsai  and  H.  T.  Hahn, 
published  by  Technomic  Publishing  Company,  Westport,  CT,  in  July  1980. 
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TAPE  #l 

PROPERTIES  OF  UNIDIRECTIONAL  COMPOSITES 


T300/5208 


Eng'g  Const. 


(A) 


Eq.  1.12 


SI -Eng 


USER  INSTRUCTIONS 

TAPE  #1:  PROPERTIES  OF  UNIDIRECTIONAL  COMPOSITES 


STEP  PROCEDURE  ENTER  PRESS 

1  For  T  300 /S20H  plies,  initialize  values 

«>f  R. ,  , Y,Y,S,F*y  ”  A' 

( =  "  1  *  >  >  ‘  y »  and  hQ  in  SI . 


DISPLAY 

1.00 


2 


Calculate  invariants  +  )  , 

.  ‘‘■id  r  stress  failure  tensors 

F  and  1  ;  strain  failure  tensors  G. 

in  i 


and  G.  . . 

».) 

Convert  from  SI  to  English  (lb,  in,°F) 
units . 


U 


5 


B 


1.00 


1.00 


3  Convert  from  English  to  SI  units. 

4  Convert  from  SI  to  kgf  (kgf,  mm,°C)  units. 

5  Convert  from  kgf  to  SI  units 

Alternatives 

0A  Clear  memory 

1A  To  initialize  other  values  of  E^,  E^ 

....  h^,  the  values  may  be  stored 
manually  in  the  data  reqisters. 

2A  Calculate  Q  ,  U's,F^,F^j,  etc. 


C 

D 

E 


CMs 

STO  18 


A 


1.00 

1.00 

1.00 


E 


x 


1.00 


Then  Steps  2-5  can  be  performed  as  appropriate. 

For  example,  one  can  initialize  in  Eng.  units  and 
convert  to  SI  by  using  Step  3. 


Computed  ply  data  should  be  recorded  in  blocks  3  and  4  for 
future  use.  Tape  #1  need  not  be  run,  unless  a  change  in  unit 
(e.g.  from  SI  to  Eng)  or  change  in  properties  is  desired. 
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Tape  #1 


Title  PROPERTIES  OF  UNIDIRECTIONAL  COMPOSITES 
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Initialize 

T300/5208 


B  SI  to  C English 
English  to  SI 


Q  SI  to 
kgf/mm^ 
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USER  INSTRUCTIONS 

TAPE  § 2 :  OFF-AXIS  PROPERTIES  OF  UNIDIRECTIONAL  COMPOSITES 


STEP 

PROCEDURE 

ENTER 

PRESS 

DISPLAY 

0 

Ply  data  must  be  in  storage 

1 

Calculate  modulus  Q. .  and  compli- 

U 

ance  S..  in  rotated  coordinate  sys- 
U 

tern  at  angle  6  (positive  counter¬ 
clockwise)  to  reference  coordinates 

8 

A 

1.00 

2 

Input  applied  stresses  in  reference 

*1 

ff2 

B 

*1 

°2 

coord,  system.  If  unit  stresses,  such 
as  [1,0,0],  are  entered,  the  resulting 

R/S 

* 

H  values  are  the  allowable  strengths. 

CT6 

R/S 

°v 

3 

Calculate  corresponding  strains  in 

-- 

C 

R 

reference  and  material  coord,  sys-  , 

1 

R/S 

r' 

terns  and  calculate  strength  ratios  R  &R 
(defined  as  the  ratios  by  which  the 
applied  loading  must  be  multiplied 
to  reach  the  failure  surface). 

Alternative  A 

2A 

Input  applied  strains  in  reference 

el 

e2 

STO  10 

61 

C2 

coord,  system 

STO  11 

e6 

STO  12 

e6 

3A 

Calculate  strains  in  material  coord. 

.. 

D 

R 

system,  and  calculate  strength- strain 
ratios  R  and  R  . 

“  “ 

R/S 

r' 

Alternative  B 

2B 

Input  strains  in  material  coords. 

ex 

STO  07 

ex 

(Step  0  needed,  but  Step  1  not  needed) 

ey 

STO  08 

ey 

€  S 

STO  09 

es 

3B 

» 

Calculate  strength  ratios  R  &  R  . 

_ 

E 

R 

R/S 

1 

R 

Steps  0,  1  and  2  must  be  executed  at  least  once 

before  Step 

3. 

If  only  the  angle  in  Step  1  is  changed  while  the 

stress  remains  the  same 

Step 

2  can  be  skipped.  If  the  stress  is  changed  while  the  angle  remains 

constant.  Step  1  can  be  omitted. 
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USER  INSTRUCTIONS 

TAPE  #3:  IN-PLANE  STIFFNESS  OF  SYMMETRIC  LAMINATES 


STEP 

PROCEDURE 

ENTER 

PRESS 

DISPLAY 

0 

Enter  ply  data . 

-- 

1 

Enter  ply  angle  9.,  and  number  of 
plies  at  that  angle,  n.,i=l,4,  and 

°1 

nl 

A 

R/S 

°1 

2  n  j 

cle»U.eV,*A.V*A,V*A,„<lV*A 

02 

B 

92 

for  each  ply  and  sum  them,  and 

n2 

R/S 

2(nj  +«2) 

calculate  n  .  Since  l<i<4,  values 

93 

C 

03 

for  nonexisting  plies  need  not  be 

n3 

R/S 

2(  nj+n2+n3) 

entered,  i.e.  for  a  [0/90]  laminate, 

94 

D 

64 

i=3  and  4  can  be  skipped. 

n4 

R/S 

n 

2 

Calculate  V..,  h,  A../h,  and  a..h. 

lA  ij  ij 

-- 

E 

1.00 

3 

Calculate  engineering  constants 

-- 

E' 

-- 

R/S 

E| 

-- 

R/S 

UI2 

-- 

R/S 

A 

Alternative  A 

1A 

Rotate  entire  laminate  by  y\  ** 

y 

A' 

2ni 

n.  remain  the  same 

i 

-- 

B' 

2(  n j  +  n2) 

-- 

C1 

2 (nj  +n2+n. 

-- 

D’ 

n 

2A 

Calculate  transformed  A./h  and  a..h. 

U  u 

-- 

E 

1.00 

3A 

Calculate  engineering  constants 

-- 

E' 

Eo 

etc. 

*  The  number  of  plies  ni  of  each  ply  orientation  are  those  in  the  upper 
half  of  the  laminate.  The  total  number  for  each  orientation  is  2nj . 

The  thickness  h  in  Register  26  is  the  total  thickness  of  the  laminate; 
the  number  in  Register  46  is  one  half  of  the  total  ply  number.  For 
symmetric  laminates,  only  the  fraction  of  each  orientation  rather  than 
the  absolute  number  of  plies  is  important.  If  more  than  four  ply 
orientations  are  needed,  additional  6  and  n  can  be  entered  through  Key  B, 
C  or  D;  but  not  Key  A  which  initializes  the  program. 

**  This  is  equivalent  to  rotating  the  reference  coordinates  in  the  clockwise 
or  negative  direction. 
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Tape  #3  I  n-Plane  Stiffness 
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345 

43 
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49  PRD 
49  PRD 


400 

15 

15 

401 

44  SUM 

402 

17 

17 

403 

44  SUM 

404 

19  19 

405 

00 

0 

4  0  6 

42  STD 
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46 
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409 
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47 

410 

4£ 

STD 
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412 
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414 
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14 
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STO 
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420 

44  SUM 

421 

46 

46 

422 

43  RCL 

423 

13 

13 

424 

^1  SEP 

425 

33  X2 

76 

TEE 

e4+r  45i 

19 

B 

10 

E  ■ 

*  452 

43 

R 

43 

RCL 

453 

20 

33 

33 

454 

42 

s 

35 

i--:=: 

455 

40 

91 

R  "S 

456 

44 

s 

43 

RCL 

457 

46 

34 

34 

4  5- S’ 

43 

F 

35 

1  -  X 

459 

19 

91 

P.'S 

460 

71 

s 

43 

RCL 

46 1 

3  3 

35 

35 

462 

00 

ercr 

- 

463 

00 

43 

RCL 

464 

0 1 j 

3  3 

3  3 

465 

00 

95 

= 

466 

QD 

94 

♦  ~ 

467 

00 

31 

PS 

46$ 

00 

43 

RCL 

469 

00 

36 

3  6 

470 

00 

35 

1  X 

471 

00 

91 

PS 

4  72 

00 

7  6 

LBL  " 

473 

00 

16 

ft* 

474 

00 

42 

STD 

475 

00 

00 

00 

476 

00 

44 

SUM 

477 

uu 

1  3 

13 

478 

00 

44 

SUM 

4  79 

00 
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USER  INSTRUCTIONS 

tape  »4:  IN-PLANE  NCNMECHANICAL  STRAINS  OF  SYMMETRIC  LAMINATES 


STEP 

0 


ENTER 

PRESS 

AT 

E 

c 

R/S 

A 

DISPLAY 

AT 
1.00 
R , 


PROCEDURE 

3( 

Data  from  Tape  #3  mast  be  in  storage. 


Enter  AT,  and  c  and  calculate 


Calculate  nonmechanical  strains  and 

the  strength  ratios  R^  for  each  layer, 
i=  1  -  4 . 

The  strengths  are  based  on  the  c ombination  -- 
of  environmentally  induced  strains 
from  both  temperature  and  moisture. 


*  Data  in  Registers  47  and  49  are  needed  for  this  Tape. 

**  This  step  is  not  necessary  for  the  strength  calculations  in  subsequent  tapes. 
In  fact,  only  side  1  of  this  tape  (program  steps  000  -  166)  need  to  be  entered. 
The  nonmechanical  strength  ratio  indicates  the  temperature  and/or  moisture 
level  that  would  induce  ply  failures. 
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USER  INSTRUCTIONS 

TAPE  #5:  IN-PLANE  STRENGTH  OF  SYMMETRIC  LAMINATES 


STEP 

PROCEDURE 

ENTER 

PRESS 

DISPLAY 

0 

Results  from  Tape  #3,  or  #3  and 
#4  must  be  in  storage* 

-- 

1 

Enter  N  components  and  calculate 

M  „ 

N  /h 

E 

N  /h 

e  .  Unit  stresses  such  as  [1,0,0] 
car.  be  entered  here.  The  resulting 

Vh 
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Vh 

S  values  are  the  allowables. 

Vh 

R/S 

1.00 

2 

Calculate  strength  ratios  R  &  r' for 

ii 
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R! 

R2 

each  layer,  i=l-4.  These  strengths 
are  based  on  the  ratios  of  mechanical 

-- 

B 

strains;  not  the  total  strains. 

-- 

C 

D 

R 

3 

R 

4 

*Tape  #5  can  be  used  following  Tape  #3  if  nonmechanical  strains  are  neglected. 
Ply  data  tape  used  in  Step  0  of  Tape  #3  automatically  make  e  =e  =0.  Tape  #4 
need  not  be  run  for  making  AT=  c  =0 .  X  ^ 


JfSf> 


TAPE  *6 

FLEXURAL  RIGIDITY  OF  SYMMETRIC  SANDWICH  PLATES 


USER  INSTRUCTIONS 

TAPE  #6:  FLEXURAL  RIGIDITY  OF  SYMMETRIC  SANDWICH  PLATES 


STEP  PROCEDURE  ENTER 

0  Results  from  Tapes  #4  for  in-plane  non¬ 

mechanical  strains  e^  shall  be  in  place. 

1 

Tape  #6  can  be  used  by  manually  in¬ 
putting  (0.,n?),  instead  of  Tape  #3. 

1  Enter  core  half  thickness  z  and  z 

c  c 

initialize  z. 

i 

2  Ply  angle  6^*and  position  z^ 

for  computation  of  V^,  V^,  V^, 

and  V  *  are  initialized.  The  thickness 
4D 

of  the  laminate  up  to  that  point  is  displayed?- 

3 

3  Calculate  V._,  12D../h  andd... 

iD  ij  U 

4  Calculate  engineering  constants 


PRESS  DISPLAY 
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2z2 
2  z  „ 

2  z , 


1.00 


*  Consistent  with  Tape  #3,  n.  refers  to  the  plies  of  0.  in  the  upper  half  of  the  laminate. 

**  Ply  orientations  are  starting  from  the  mid  plane,  z  =  0,  or  the  upper  surface  of 
the  sandwich  core,  z  =  zc;  this  is  the  opposite  of  the  usual  ply  orientation  code 
which  uses  an  ascending  order  from  the  bottom  ply,  z=  -h/2. 
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Title 
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USER  INSTRUCTION 

TAPE  #7:  FLEXURAL  STRENGTH  OF  SYMMETRIC  SANDWICH  PLATES 
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Results  from  Tape #6  must  be  in 
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Calculate  ply  strains  and  the  strength 
ratios  at  the  outer  face  of  each  ply. 
These  strength- strain  ratios  are  based 
onmechanical  strains,  not  the  total 
strains. 
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PLY DATA 
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Ply  Data  Tape 


USER  INSTRUCTION 


TAPE  #8:  THE  V'S  AND  NONMECHANICAL  FORCES  OF  GENERAL  LAMINATES 


STEP 

PROCEDURE 

ENTER 

PRESS 

DISPLAY 

0 

Ply  data  must  be  in  storage 

• 

1 

Enter  total  number  plies. 

n 

A 

0 

2 

Enter  ply  angle  0^,  and 
number  of  plies  at  that 

61 

R/S 

01 

angle  n^ . 

nl 

R/S 

V«1 

6  2 

B 

,  n  ,  . 

(T+nl> 

n2 

R/S 

V<°2 

03 

C 

(-?+nl+n2> 

n3 

R/S 

V«3 

94 

D 

|-7+nl+n2+n3> 

n4 

R/S 

V«4 

!5  • 

B** 

(-J+n1+. . .) 

* 

* 

* 

3 

3 

Calculate 

E 

4/n 

4 

Calculate  N?*,  M^* 

AT 

i 

E 

AT 

c 

R/S 

mN* 

M6 

** Additional  plies  can  be  entered 

through  B , 

C,  or  D 

key. 

Key  A  cannot  be  used  for  this  purpose  because  it  initializes 
the  summation  process. 
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USER  INSTRUCTIONS 


TAPE  #9:  MODULUS  OF  GENERAL  LAMINATES 


STEP 

0 

1 

2* 


PROCEDURE  ENTER 

Results  of  Tape  #8  are  in  storage. 

*  *  * 

Calculate  A. B. D. . 

i]  13  13 

Calculate  A^ ,  B^,  D± ^ 


PRESS  DISPLAY 


26 


B 


26 


*  This  is  a  necessary  step  before  the  matrix  inversion  in  Tape  #10. 
Inversion  of  the  normalized  modulus  in  Step  1  is  not  meaningful. 


49 
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